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A multiple-time scaling analysis of the dissipative, transversely driven Landau-Lifshitz equation in 
presence of exchange, shape demagnetisation and week anisotropy fields is performed for a dynamic 
domain state. Stationary solutions of the resulting equations explain the spatiotemporal structure 
of the walls and are in agreement with previous simulations. 
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I. INTRODUCTION 

in : 

' In the interdisciplinary, unifying work on structure formation in driven dissipative systems ferromagnets are 
only treated in exceptional cases. This is surprising as ferromagnets exhibit already in static fields domain structures 
which are well understood (for reviews see 1^) and which even serve as textbook examples Dynamical questions 
have also been investigated. This has been done by extending the research [|j to weak driving fields (compare e.g. 
.^fH ■ Q). For stronger driving fields numerous work (reviewed in [^) has been published to describe the ferromagnetic 
resonance instabilities 

Apart from single exceptions like 0, where in simulations dynamic domains for a special model were reported, the 
, existing investigations on dynamical questions are limited to the linear or to the weak nonlinear regime. An approach 
I ' 1^ to analyse the strongly nonlinear dynamics of driven ferromagnets in terms of structure formation has recently 
J> , been proposed by the present author. In this approach a transversely driven model was investigated containing an 
' exchange, an anisotropy and a shape demagnetisation field. Numerical simulations exhibited in the rotating frame 
, a stationary domain structure with a precessing motion in the wall regimes. Some of the characteristic elements of 
' this structure have been explained in Q by analytical methods. The spatiotemporal dependence of the domain walls, 
however, remained unexplained. 

It is the aim of the present investigations to work out analytically an explanation of these dynamic domain walls. In 
0^ , sec.|^ the basic ingredients of the present approach are given. After performing a transformation in sec. Ill a multiple 
time-scaling analysis is employed in sec.IV to deduce reduced equations of motion on the slow time scale. Representing 
^ , the general result of this work these eqs. are in sec.0 applied to determine the spatiotempo ral d ependence of the 
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stationary walls. In sec.0 the results are compared with the simulations of Finally in sec. VII 
are given. 



some conclusions 



II. THE MODEL AND BASIC ASSUMPTIONS 



X 

5_] At a mesoscopic scale the dynamics of ferromagnets is governed by the Landau-Lifshitz equation, which takes the 

form 

9t m = —(1 + Fmx) mx {h^^ez + h±ex — m + J Am -\- Anizez) + ujmxez (1) 

in the frame rotating with the driving frequency uj around the direction. m(r, t) is the local magnetisation in 
the rotating frame being related to the magnetisation in the laboratory frame miab by miab = exjp^utezX ) m. 
The amplitudes of the external static and the external circular driving rf field are denoted by and by h± , 
respectively. The term m ~ J mdV represents, in reduced units, the demagnetisation field of a sphere of 
volume V. The contribution JAm results from the isotropic ferromagnetic exchange interaction and the uniaxial 
anisotropy is described by ^m^e^ where ^ < is assumed in this work. The Landau-Lifshitz damping rate is 
represented by F. The gyromagnetic ratio and the magnitude of the magnetisation m — |m| are set equal to 1. 

A is assumed to be a small quantity which permits an perturbative treatment. Recall that eq.(^) results from a 
continuum approximation in space (compare Q) which implies a slow spatial variation of m. Thus the exchange field 
JAm is also treated as a perturbation. Introducing the expansion parameter e = —A > and scaling the position as 
r (—J/Ay^'^r the perturbational fields of eq.(^ JAm -I- ^m^e^ take the form e(Am — m^e^). 

Due to the rn{t) term, the problem described by eq.(^ is of the mean field type and the usual technique can be 
applied. This is to employ an ansatz for m, solve formally the problem for given m and in the last step examine the 
results for self-consistency. 
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Based of the findings of |^ as ansatz a stationary domain structure is assumed in which the magnetisation is 
constant nearly everywhere and takes only two values m_|_ or m_ realised in the generally disconnected partial 
volumes V+ — n+V and V- — n^V , respectively. The wall regimes where the magnetisation changes from m_ to 
m_|_ are assumed to be narrow. Then the volume of the wall regions can be neglected compared to the sample volume 
V which implies 

nidom = n-|_m_|- + n ; with n+ + n_ = 1 . (2) 

This work is focusing on the stationary domain solutions of eq. (|I]) . To keep the analysis as simple as possible m(t) 
is replaced by nidom from the beginning. For a justification of this ansatz it is pointed out that already for e = 
relaxation of m(t) to mdom was found in |^ . 

In this previous work the zeroth order contributions to n± and to in± have been calculated. The perturbation 
theory of this work makes it necessary to separate in notation between full quantities and their contributions to the 
e expansion. This is archived by generally setting a = a + ea*-^-* + O(e^) for an arbitrary quantity a. This rather 
unconventional notation of the zeroth order contributions avoids complicated indices. Employing this notation the 
results of Q take the form 

- ± ; m± = It (Twei + Th^ea) (3) 
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d = ^-u; ■ u=(d2 + r2)i/2 . v = {u^-T^hlf^^ ; w ^ {u^ + hl)'/^ (4) 
and where the transformation from the internal orthonormal system (ei, e2, 63) to (e^;, ej^, e^) is found to be given by 

(5) 

With eqs.(|3|-|^) and with the above replacements eq.(|l|) takes the form (compare [^) 

dtm = —ujw^^ mx ( uei + Thj_nixe2) — e (1 + rmx)mx (Am — m^e^ — m|^^^j^) + O(e^) (6) 

where m^^j^ denotes the first order contributions of mdom according to the convention in notation. 

Eq.(||) is the starting point for the further analysis. Compared to eq.(|l|) this eq.(|^) is constrained by m(t) = nidom- 
Note, however, that eq.(|^) still describes in general nonlinear deviations from the domain state as long as these 
deviations are local. 

The m± have to satisfy eq. (||) with 9t m = Am = from which 

m« = m±x[r - m±x] [±-m^l - —e, + ^e^] (7) 

UIV UOW UJW 

results. Eq.(||) leads in first order to 

— (1) (1) / - ~ \ , - (1) , ~ (1) ro\ 

^dom — (™+ ^ rn ) + n+m]^' + n_m_' (8) 

where the ih± and the n± are given by eq.(^). Note that from eqs.(^ and (||) the dependence of m^^j^^ , m^'' and of 
n^'' on the model parameters h\\, h±,uj and T is not completely determined. For such a determination an additional 
relation between these quantities is needed. Thus at this stage of the calculation m^^^^j^^ and m^'' are treated as 
fimctions of n'^P and of the model parameters. 



III. TRANSFORMATIONS 

In zeroth order the dynamics described by eq.(j^) is hamiltonian Thus the dynamics should be formulated in 
variables, which in the e = limit reduce to angle - action variables. The total transformation to these variables 
is divided into a sequence of single transformations. First an orthonormal system (61,82,63) is defined by ei = 
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(2 6i) -"^(ni- — ni+) by 63 = (263) ^(m_ + m+) and by 62 = 63X61 with bi = — m+m_)/2 and with 63 = 

a/ ( 1 + m_|_m_)/2 . As the m± are the full domain magnetisations, the definition of the (61,62,63) is independent 
of the perturbative treatment and holds to all orders in e. Note that are the zeroth order contributions to 6^ in 
accord with the general convention. 

In the next steps this (61, 62, 63) base is first rotated about the (2)-axes and thereafter standard spherical coordinates 
0, $ with respect to the new (3)-axis are introduced. These transformations are explicitly described by 

sin Q cos $ sin 8 sin $ cos 6 \ / cos 6 

cos 8 cos $ cos 8 sin $- sin 9 1 | | ^2 | (9) 

-sin$ cos$ / I 63^(6) 





where p(8) is given by p(8) = yb^ — cos^ 8 . Finally a change of variables from 8, $ to new independent variables 
Q,(p is performed by setting 

tan ^ ^' = — — -tanf7- with = 1/63^ - 1. (10) 

2 p(8) - sin 8 2 V 3 ^ ^ 

which completes the total transformation. 

Eq.(^ implies d$ = p"^(8) cos8sin$d8 + .g(8,(^)diy9 where the function g(8, (^) is defined as (7(8,1/?) = 
p(8) + sin 8 cos With this expression for d$ eq.(||) leads to 9t m = g p^^ dt Q + g sin 8 6$ (9f . To transform 
the rhs. of eq.(|) in a first step ei and §2 are expressed by the 6^ using eqs.(^,(|^) and the definition of the 6^. 
After that eq.(p) is employed to represent the rhs. of eq.(^) in terms of 69 and 6$. Introducing a scaled time i by 
t = Lo u^bib^v t and setting e — rh±ujw~^e the results of this basically straightforward calculation leads to 

di 8 = ep(8) {Xl^ + X|" + rxg^ + rxg" + Ye) + O(e') 
sin 8 {df ^-l)= i (-Xg' - X§' + TXl^ + rX|" + Y^) + 0{i^). (11) 

The quantities g Xq^ = 6eAm and g X^ = 6$ Am are calculated to 

^Xg' = V(p"\9V8) - [p-isin8V8 + g COS8V93] (gsin8V^) 

gXI'' = V(5sin8V(/>) + [p^^ sin 8 V8 + .9 cos 8 V^j] (p-^g V8 ) . (12) 
From the definitions g Xq^ — — (m62) (696^) and g X'^ — — (m62) (6^6^) 

gX^ = /i^w"^(-sin8sin$ + r5 cos8)[cos8sin<I> + r(5sin8 - cos$)] 

gX^ = h\w^'^{- sinQsin^ + Vg cos 8)(cos $ + F cos 8 sin $) (13) 

is obtained and for Yq and 1$ 

Ye ^ -T^h^v-^[u~em^2^^, + h^-e^m^Z^ sin 8 

-F/i^w"2^"i[(l-F2)cose(pcos'I> +sin8) + 2rpsin$] 

^Th^v-^u-\{v'' -V'^hDeim^ll^ + V'^h^ue^m^^^J sin8 

+ F/iiw"2g"^[-2Fcos8(pcos$ +sin8) + (1-F2)psin$] (14) 



results with $ $(8, (/>) as defined by eq.(|lO|). The quantities p and g stand for the zeroth order contributions to p 
and to g and are explicitly calculated to 



p(8) = V^^Tsii^e , g(8,(p) =p(8) +sin8cos$ with n^-^. (15) 

F/i I 



In the zeroth order limit eqs.(|ll[) reduce to 9j 8 = and to d^ip = 1 . Thus as required the variables 8, ip become 
angle - action variables in the 6 = limit. A second feature of the transformation should be pointed out. According 
to the definition of the 6; and to eq.(||) m(8 ^ 0,ip) — m. and m(8 — > tt, (/?) = m+ holds. Thus the values 8 = 
and 8 = TT characterise the interior of the domain regimes and a domain wall corresponds to a change from 8 = to 
8 = TT. Note that, as it must be, eqs.(0) are satisfied for 8 = and for 8 = tt with df 8 = V8 = 0. 
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IV. MULTIPLE-TIME SCALING 



Next the standard multiple-time scaling method is applied to the eqs.(|l l|) . Denoting the slow time scale asT = it 
the ansatz 

e(£,r) = e{T,r) + ie^^Hi,T,r) + 0(e2) 

ifii v) = i+ ^{T, r) + e>(i) (t, T, r) + 0{e^) (16) 

is employed which already satisfies the zeroth order of the perturbative expansion. In first order the t-secular terms 
of e(i) and of (ys^^) are as usual eliminated. This procedure leads to 

drO^ m ( (Xl^) + + r(Xg^) + T{Xg^) + (Ye) ) 

sine dr^^- - (xg") + Tixi'^) + r(x|") + (y$) (17) 

which determines the T- dependence of e{T, r) and of if{T, r). In eqs.® ( /(O, ip) ) = fi/(27r) /g^''" ' /(O, < + (^) dt 
denotes an averaging of an arbitrary function f{Q,(p) over the fast time variable t. Using a substitution analogous to 
eq.(|l^) these averages can be performed and result in 

(Xg^) = (p + n)V —^-^ 17 sin e cos 6 (V ipf 

ra = sme^v(^v^) 

p — \l p + il 

{xg') ^ hlw-^ {n - p) [1 ~r^{i +np + n'^)] cote 
= 2 r hlw-^ {n - p) cot e 

{Ye) =-^2;^^«-2[uelm|,l + /^xeam^l] sin 9 

(Kf) =Th^v-^u-^[iv^ -T^hDemi^l^ + T^h^ue^m^^lj sin 9. (18) 

where now p stands for p(9). 

Together with the expressions ( |l4| ) the coupled set of eqs.(p7|) represents the general result of this work. On the 
slow time scale T the temporal and the spatial evolution of the magnetisation is governed by these reduced equations 
of motion. It is expected that various, interesting questions can be analysed on the basis of these equations. One of 
these problems, the profile of a stationary domain wall is treated in the next section. 



V. DOMAIN WALL 



In this first approach to the domain wall problem I restrict myself to planar walls whose spatial dependence varies 
- locally and/or approximately - only in one direction, the normal direction of the wall. The analysis is further 
restricted to stationary walls. Thus 9(^) is assumed to be time independent. For the cyclic variable (p the ansatz 
<^(T, ^) = O'^' T + (^(^) is used where fl'--^"^ describes a first order frequency contribution independent of ^. Fixing 
the scale of ^ by the replacement (re^) — > h±w^^£^ and introducing the ' local wave number variable ' by 

Q(e)-^'(r,o (19) 

eqs.([l7|) lead to a system of ordinary differential eqs. for 9(^) and for Q{£,) 

{p ~ T^{1 +Qp + Cl^)] cot 9 + Ce sin 9 

2r(p- fi) cot9 + C.E, sin9 (20) 



where p still stands for j5(9) = Vil^ + sin^ 9 and where the constants Ce and C$ satisfy the linear system 



(p + n) 



9' 



p{p + n) 



— tl sin 9 cos 9 Q — 



. ~ p + ri f p — 
sm 9 ■ 



» — n VP + ^ 
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rCe + C$ = r\;2/,-i^-2[^e^;fjj«^ + h^e^mll^^] (21) 
-Ce + rC* = w''hl^ - rw''hl\-^u~^[{v^ - r^/i^) eim[,'j„^ + r^/i^ueamW^J. (22) 

Apart from solving eqs.(|2^) O(^) and Q{^) have to satisfy at the boundaries 

e(C^-cx)) = 0; e(C->oo) = 7r; e'(C -> ±cx)) ; Q ^ ±oo) ^ const. (23) 

where the range of ^ has been extended over the entire interval from — oo to oo. Just by counting constants of 
integration it is obvious that the conditions (^) can generally not be satisfied, noting additionally that one constant 
is the center of the wall due to translation invariance of eqs.(|2^). For special values of Cq and C*, however, 
solutions of this boundary eigenvalue problem may exist. 

For Ce = C$ = a solution exists which exhibits the symmetry properties 6(0 = tt — 6(— ^-^d. Q{£,) — Q(— C)- 
This will be demonstrated by construction and for this propose eqs.(pO|) are formally integrated 




Q=2T^--^ dp , ' . (24) 

Consider both & and Q as functions of O or according to the relation sin 6 ~ (57^ — p^)^/^ as functions of p. Then 
the eqs.(p^ represent a system of integral eqs. for & and Q in dependence of p. Denoting the explicit solutions of 

this system by F{p) = Q' and by G{p) — Q further integration leads to ^ = j^^2 P~^W + sin^ ©) d6 where the 

center of the wall has been fixed by 8(0) = 7r/2. Thus the ^ dependence of is implicitly determined from which 
Q{C) and consequently (p{T,S) can be calculated according to eq.(|l9|). 

In the limit p ^ which corresponds to the limit ^ ±oo the eqs.(p4|) lead to F{p) ^ {p — f^)^/^ and to 
G{p) —^ const. . These findings are valid for all parameter values of 17 and of F and imply that all the boundary 



conditions (23) are satisfied. 

Turning now to the question whether further solutions of the boundary eigenvalue problem exist eqs.(^0|) are again 
integrated formally. Due to Ce ^ and C$ ^ additional contributions arise to eqs.(p4|). From these contributions 
the asymptotic behaviour for ^ -> oo (e'(-O)^ - {Q'iOf ^ Ce and Q{-^) - Q{^) - C$ (p - Cl)^'^/^ can be 
deduced. These results imply that the conditions ( |23| ) can not be satisfied with Ce 7^ or with C$ ^ 0. Thus there 
are no further solutions and Cq = C$ = represents a compatibility condition. 

In general the explicit solutions of eqs.(pO|) with Cq — C$ — have to be calculated numerically. Fig.(l) shows 
the results of such a calculation for the parameter values tl — and F = 0.1 . Note that the values of Q are small 
compared to the variation of 9. This findings imply that the spatial dependence of the wall is mainly determined by 
&{£,) and are attributed to the fact that Q{^) is proportional to F. Thus alternatively an iteration procedure can be 
employed to find approximate solutions. Starting with Qo{0 = at the rhs. of eqs.(p^) these eqs. determine first 
approximations. For the above parameter values the deviations from the exact solutions are found to be small. As 
usually F <C 1 is satisfied this indicates that the first iteration of eqs.(p^) - which can be calculated analytically - may 
be sufficient in many cases. 

According to eq.(pT|) Cq = C$ — implies the relation ueim"^^^ + h± e2m^^^^^ ~ 0. This additional relation to the 
6qs.(0) and (^) now determines in a unique way the dependence of n^'' of m"^^ and of m^^jj^ on the model parameter. 
With eq.(p2|) this also applies to il'^'. Thus all quantities of the domain states including the spatiotemporal wall 
structure are uniquely determined. This result together with the fact that the wall width remains finite in the large 
V limit shows the self-consistency of the ansatz(||). 

A further consequence of the additional relation is (Yq) = and thus this term in eqs.(p^) drops out. By a 
redefinition of the fast time scale in the first order term, it is possible to transfer the contribution (1$) of eqs.([l7|) 
to the fast dynamics. Employing this redefinition both (Yq) as well as (y$) drop out of eqs.(p^). This findings 
considerably simplify the analysis of further problems based on the general eqs.(p7|). 

Next the internal fields Hi and H2 are introduced by H2 — Hi +Luez = h^^e^ + /i±e^ — rndom- Employing again 
the above relation and eqs.(||-^ it is found that up to the first order 

HiH2 = (25) 
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holds which generahses the zeroth order result of 1^. In this ref. the importance of eq.(p5|) as a criterium for the 
existence of the domain states has already been discussed. 

To conclude this section it is mentioned that within the domain regions a linear stability analysis has been performed 
with the result that all the in zeroth order undamped modes of have now finite damping constants proportional to 
e. This stability analysis does not include the wall regimes as for such a complete treatment numerical methods have 
to be applied. The complete zeroth order stability analysis, the present partial first order analysis and the results of 
the computer simulations yield considerable evidence for the stability of the dynamic domain states. In this context 
it is pointed out that in the weak coupling limit F <C 1 a Lyapunov function can be constructed on the slow time 
scale which implies stability. 

VI. COMPARISON WITH SIMULATIONS 

With the numerical results of the last section and an additional integration the zeroth order contributions of eqs.(|l6|) 
are determined. With these results for Q and if the Cartesian components of rui can be calculated according to eq.(|^. 
This has been done for the parameter values = = 2, hj_ = 0.5 and T = 0.1 neglecting H,^^") and all other first 
order corrections in eq.(^) which implies also = e^. The results of this calculation are given in Fig. (2a) and Fig. (2b) 
where the spatiotemporal dependence of cos0 = mi(m^ +7713)^^^^ and of 7713 are presented. For the time dependence 
anharmonic oscillations are found which are plotted on the t — ni scale for one period. 

Fig. (2c) and Fig. (2d) show the corresponding results (D of the computer simulations of eq.(||) performed with the 
same parameter values and using in addition J = 0.01 and A — —0.005. The original simulation results are rescaled 
to the theoretical scales. With the exception of the center of the wall and the time origin no fitting procedure has 
been used. 

With deviations of usually less than one percent the analytic results agree with the simulations. Quantitative 
agreement for rildom, for m± and for was already found in ||^ with a similar high accuracy. Based on these 
agreements it has been demonstrated that the employed methods are a suitable tool to explain the characteristic 
features of the dynamic domain states. 

VII. CONCLUSIONS 

In this work a rather realistic model for a ferromagnet was investigated. Depending on the parameter values this 
model can exhibit stationary dynamic domain states in the rotating frame. An important element of this states is the 
spatiotemporal structure of the domain walls. With the present investigation this structure is understood to a high 
extent. 

The analysis of wall structure is based on the reduced equations of motion (0) which from the theoretical point of 
view represent the main result of this work. These eqs. which govern the temporal and the spatial evolution on the 
slow time scale should be applicable to other phenomena. Obvious extensions are investigations on drifting walls, on 
multi-dimensional walls and on the interaction between domain walls. 

As a further propose the present approach should be extended to take stray fields into account. This would imply 
that the dipole interaction is completely included and the model will become very realistic for ferromagnetic materials. 
Moreover by analogy with the static case, it is expected that these dipolar stray fields contribute crucially to selection 
mechanism for the wall positions. In this context the formation of regular patterns in driven ferromagnets seems to 
be possible. 

The author has benefited from discussions with W. Just and G. Sauermann. This work was performed within SFB 
185. 
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FIG. 1. Numerical solution of eq.(20) (with Ce = C"* = ) for ©(^) and Q{^) satisfying the boundary conditions, of 
eq.(23). The parameter values axe Cl = 3^^^ and F = .1. 

FIG. 2. Spatiotemporal dependence of a stationary domain wall in reduced units ^ and t. In Fig. (2a) and (2b) the results 
of the analytical treatment of this work are presented. For comparison Fig. (2c) and (2d) show the results of the numerical 
simulations of [8]. For both cases co8(j> = mi(mf + m'i)^^^'^ and ms are plotted, where rrii are the components of the local 
magnetisation in the internal coordinates, defined by eqs.(4) and (5). The scaled position ^ and scaled time r are given by 
^ = ■whJ^{—A/J)^^'^ e^r and hy t = Clt = ujvw~^ t, respectively. Compare text for the parameter values. 
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